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Abstract 

Wc identify the scaling region of a width (^(n" 1 ) in the vicinity of the accumulation points 
t = ±1 of the real roots of a random Kac-like polynomial of large degree n. We argue that 
the density of the real roots in this region tends to a universal form shared by all polynomials 
with independent, identically distributed coefficients Cj, as long as the second moment a = E(cf) 
is finite. In particular, we reveal a gradual (in contrast to the previously reported abrupt) and 
quite nontrivial suppression of the number of real roots for coefficients with a nonzero mean value 
fj, n = E(ci) scaled as [i n ~ rT 1 ! 2 . 



1 Introduction 



Recently, there was an essential resurrection of interest in statistical properties of zeros of random 
polynomials and random analytic functions [T]-|15j. 

The distribution of real roots for random polynomials was studied in classical papers by Bloch and 
Polya [TJI|, Littlewood and Offord ^7j> Kac ^Hl an d many others. Many results and an overview of 
the field is given in the recent monograph by Farahmand 

As is well-known, for the case of algebraic polynomials 



fn(t) = c + cit + ... + c n _ii 



n-1 



(1.1) 



where the real coefficients Cj are independently and identically distributed around mean zero, the roots 
are mostly located near ±1. A natural intuitive picture for such an accumulation was proposed by 
Kac. Indeed, if all coefficients Cj are of the same order then only for \t\ = 1 will the terms have a good 
chance to interact and cancel each other producing zero. 

In particular, for Gaussian-distributed coefficients with mean E[cj] = and variance E[dj] = a, 
the mean density of real zeros of an algebraic polynomial p n (t) is given by the Kac formula 

y/A n (t)C n (t)-B*(t) 



Pn(t) 



nA n (t) 



(1.2) 



where 
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A n (t) = a^2t 2j = a- — -, B n (t) = -— A n (t) , C n {t) = --^A n {t) + -—A n {t) . 
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Usually, one is interested in the limit of polynomials of high degree n. For a fixed t the density is 
given by: 

(1.3) 



lim p(p(t + ±1) =p {f \t) 
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and 



Pn(±l) = - 
7T 
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The expected number of real zeros is known to have the following asymptotics: 

2 



A" n (-oo,oo) 



p n (t)dt = -\ogn + C + 0{n~ z ) 
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(1.4) 



where the leading logarithmic term was found by Kac himself, and the corrections, in particular, the 
constant 

2 



C 




1 
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sinh 2 v 



dv 



(1.5) 



was found by Wilkins [20]. In fact, the asymptotic result N n (—oo, oo) = ^ logra was proven to be very 
universal, i.e. insensitive to the details of the distribution of coefficients q, see discussion and further 
references in [22] • A similar universality result is also available for the variance of the total number of 
real zeros |21) . 

The nature of the informal arguments by Kac suggests that close to the values t = ±1 many 
random terms of comparable magnitude add up to form the polynomial, and the emerging object 
should be universal in the spirit of the central limit theorem. 
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Another natural question arises about the expected number of real zeros for coefficients with a 
nonvanishing mean value E[cj] = fM. It is straightforward to include the mean value to the stan- 
dard derivation of the Kac formula for gaussian-distributed coefficients and arrive at the following 
expression: 

K(i)=Pn(t)e-^ [El(i)+E2(t)1 y dg?cosh( g v^(t))e"^ (1-6) 
_ A n {t)Dl{t) _ Gl{t) 

D n (t) = -G n (t)-^ + -G n (t) , G B ( t ) = M £V = M .— (1.8) 

and P n (t), A n (t),5 n (t) and C n (t) are defined in (fL2j) . 

Earlier investigation [22] has revealed that any, whatever small, mean value \x > asymptotically, 
for large n, kills exactly half of real zeros and converts the above quoted result ITTil for their number into 
N n (—oo, oo) = - logra. It is this surprising abruptness that motivated us to have a closer look at the 
origin of such a behavior. This question turned out to be also related to exploring the appropriately 
scaled vicinity of the accumulation points. 



where 



2 Scaling and Universality 

It is quite clear (and can be verified numerically) that the root density function p n (f) as well as its 
limiting form 1)1.3(1 are both not universal. On other hand, the essence of the Kac argument discussed 
above and the mentioned universality of N n (—oo, oo) suggests that universal features have to emerge 
in close vicinity of the accumulation points t = ±1. In fact, it is easy to understand that the relevant 
vicinity of \t\ = 1 is of the order of n~ 1 for large degree n. We will refer to such domain as the 'local 
scaling regime', as opposed to keeping the distance |1 ± t\ from the accumulation points t = ±1 fixed 
in the limit n — > oo, the latter regime being referred to as the 'global' one. In particular, our main 
goal is to verify Proposition 12.11 

Proposition 2.1. Let coefficients {cj}J~ of a random polynomial f n (t) be i.i.d real random variables 
with E[cj] = fi n > and finite variance chosen to be unity: E[c|] = 1. Let N n (xi,X2) be the number 
of real zeros of f n {t) in the interval [l — , 1 — ^1 , with < x\ < X2 < n. Then 

lim E [N n (xi, X2)] = / p a (v)dv (2.9) 

n ^°° J XI 

where a = lim n _ >00 (n/i^) and the density function p a (v) is given by 



1 roo 2 

Pa(v) = -\H--^ e ~ aJ(v) dqq cosh (q^/c7M^))e- 3 2 (2.10) 
7r V v 2 smh v Jo v / 



sinh (v/2) 1 . . sinhv/w 1 — s'mhv/v 

{v/2) 2 1 + sinhw/w ~ cosh 2 (v/2) 1 + sinhu/v 



A similar statement is valid for the vicinity of second accumulation point t = — 1 for fj, n < 0. 
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Here we present explicit arguments in favour of the validity of such a proposition for the simplest 
case u n = 0. The 'local regime' formula ()2.1(Jj) reduces in that case to the following simple expression: 



po(v) 



1 / 1 



1 



2n V v 2 sinh 2 v 



(2.12) 



The modifications required to include nonzero mean a n are self-evident and left to the reader. 

Our starting formula is the following representation for the number of real roots N n (xx, x 2 ) of f n (t) 
in the interval [l - a , 1 - £l] , see e.g. P 



where 



N n (x 1 ,x 2 ) 



T> n (x,y;t) = E 



1 _ ZL 



1_ ±2l 



dt / dy \y\T> n (0,y;t) 



n-l 
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(2.13) 



(2.14) 



stands for the joint probability density of the random polynomial f n (t) and its derivative over t. Let 
us introduce in the above representation the scaling variables v, x and y by relations: 



t = 1 — v/n x = n 1//2 x y = n ilA y. 



,3/2,-, 



Changing the variables of integration one finds: 



N n (xi,X 2 ) 



•('2 



dv 



dy\y\V n (0,y;v) 



where 



V n (x,y;v) = E 



5 x 



— J^Cjil-v/ny 
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(2.15) 



(2.16) 



Evidently each of the scaled random variables x and y is a sum of n independent, although 
not identically distributed terms, with the magnitude of fluctuations of each term depending on the 
summation index j. Calculating the variances of the scaled variables we find that they are given by 



n-l 



A n (v)=E[x 2 ] =±J2( 1 -v/n) 2j 



and 



n-l 



B n (v)=B[f] =-L^ j (l_ t ;/n)^- 1 
and their covariance is just 



1 - (1 - v/n) 
v(2 — v/n) 



1 d (a n -\ 
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C n (u) = E [xjf] = — V j 2 (l - „/„ 
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(2.17) 



(2.18) 



(2.19) 



We further notice that all three quantities A n (v), B n (v) and C n {y) have a well-defined finite large n 
limit: 



A n (y) -> A c 



1 - e 



-2c 



2r 



1 dA. 



2 cfo 



C n (v) - C ( 



1^4« 

4 <fo 2 



(2.20) 
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coinciding with the limiting values of A,B,C in Eq. (|1.2|) after rescaling t = 1 — v/n and the limit 
n — > oo for a fixed v. Now invoking the local central limit theorem [2f)j we infer that the limiting 
joint probability density of the scaled variables x & y tends when n — > oo to the normal law 
with variances , and covariance , as if the polynomial coefficients were normal with unit 
variance. Therefore, the Kac formula for the density of real roots should be asymptotically valid in 
this case, and the formula Ea. (|2.12|) immediately follows after substituting those limiting values into 
the Kac expression Eq. (|1.2() . Separation of the two regimes and identification of the scaling is in some 
sense similar to 'local' versus 'global' scaling regimes in spectra of random matrices. Only when the 
spectral parameters are scaled appropriately will the various correlation functions characterising the 
eigenvalues of large random matrices show a surprisingly robust universality |24[ I25j . 

Although the expressions Ea. (|2.1U|) . and especially Ea. H2.12|) look very natural (cf. the structure 
of the asymptotic result Eq. (|1.4jl ) and their universality stems from such basic fact as the central limit 
theorem, we failed to trace a similar statement in the available literature on random polynomials with 
real coefficients. For complex zeros a kind of 'local' regime was studied in much detail by Shepp & 
Vanderbei |3j and Ibragimov & Zeitouni [3| who showed that those zeros tend to concentrate asymp- 
totically on the unit circle. Only very recently (squared) expression of the type l|2.12[) independently 
emerged in studies of complex zeros of polynomials with complex i.i.d. coefficients, see Shiffman and 
Zelditch, ^Sj. A different kind of universality results was addressed in a very recent paper [7j which 
appeared when the present paper was under completion. 

3 Numerical Examples, Discussion and Perspectives 

It is interesting and informative to have a look at the profiles of the real root density in the 'local scaling 
regime' around the point t = 1 (i.e v = 0) for some typical values of the parameters, and compare 
them with the results of numerical simulations. In fig. JJ) we plotted the density profiles obtained by a 
direct numerical search for real roots of polynomials with gaussian-distributed coefficients for degrees 
n = 100, 300 & 1000. The (scaled) mean value a n = fj^n of the coefficients was chosen to be a = 10 
and the variance was always kept unity. The results can be compared with both the exact predictions 
of the Kac-type formula (|1.6|) for the same values of n, a and \i n = y/a/n, and with the asymptotic 
universal profile of (|2.10j) . It is trivial to see that as n increases the exact Kac-like expression begins 
to coincide with the asymptotic profile, and the analytical curve agrees well with the numerical data. 
To verify universality we also performed numerical simulations for the case of coefficients uniformly 
distributed in the interval of the widths l/Vs around the same mean values fj, n . The picture looks 
very similar, see fig. (J5J), and again agreement with the asymptotic formula for large values of n is 
very good. 

One can notice a couple of interesting features. First simple observation is that the exact Kac 
formulae are asymmetric about v = for 'small' n (in the 'local regime') but it becomes "asymp- 
totically symmetric" with increasing n, as predicted by (|2.10j) . This feature is easy to understand 
in view of the global inversion symmetry t — > t" 1 which holds exactly for all polynomials with i.i.d. 
coefficients. In the local regime close to accumulation points this symmetry implies the asymptotic 
reflection symmetry v — > — v. 

The most surprising feature is a non-trivial double-peak structure of the density profile, and it 
deserves to be discussed in more detail. 

In fact, it turns out that as a increases the shape of the density profile described by formula (|2.10j) 
changes from that with one maximum to that with two symmetric maxima as illustrated in the fig. 
©• 
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Figure 1: Numerical vs. analytical density profiles for real roots of polynomials with gaussian- 
distributed coefficients. Each data set corresponds to 100, 000 realisations of random polynomials. 



This fact can be verified analytically by considering the small-v expansion: 
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(3.21) 



valid for av 2 <C 1. We see that for a > a c = 6/5 the local maximum at v = is converted to 
a minimum. On the other hand, for |u| > a we easily find from Eg. (|2. 10(1 the a— independent far 
tail decay law p a {v) ~ 1/V|u| which ensures that at least two symmetric maxima have to appear at 
?) / 0. Notice that the far tail values are much larger for a>l than the exponentially small density 
Pa(0) = e _Q//2 at the origin, and thus the maxima have to be quite pronounced. 

The precise reason for such a nonmonotonous behavior is not clear to us at the moment, and deserve 
further investigations. We consider that feature as an indication of a rather nontrivial statistics of 
zeros accumulating in the scaling region around the points dbl. 

Fixing the value of \i and performing the limit n — > oo amounts to letting a — > oo. This results 
in vanishing density of the real roots p a {v) — > for any fixed v. This fact corresponds to observed 
suppression of exactly half of the real zeros (the vicinity of the second accumulation point is not at 
all affected by the nonvanishing mean fi n ). In contrast, correctly tuning the mean value with the 



parameter n by letting [i Ti 



n 



~ l l 2 and magnifying the vicinity of the accumulation point via the 



correct scaling 1 — t ~ n 1 results in a density function (|2.1U|) and thus describes a gradual (as 
opposed to abrupt) suppression of real zeros in the scaling regime, see fig. (j3J). 
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Let us now briefly discuss open questions and perspectives for further research. 

A very interesting development of the theory of random polynomials comes from the paper by 
Bleher and Di 0. The authors discovered nontrivial correlations between positions of real zeros of 
algebraic polynomials. The correlations can be conveniently characterised by correlation functions of 
the polynomial zeros: 



K(t!,...,ti) = lim E 



n^oo 



N(f\t 1 ,t + A 1 )...N(f\t l ,t l + A l ) 

hm j-r— : j-r— : 

Ai,...,A,-0 Ad... A; 



(3.22) 



where N^'(a, b) is the number of real zeros of the polynomial f(t) in the interval [a, b]. 

The consideration of correlation functions studied by Bleher and Di was restricted, in our terminol- 
ogy, to the 'global regime'. It should be possible to derive 'local regime' formulae for those quantities. 
The expressions are expected to be universal in the same sense as the 'local regime' formula for ex- 
pected number of roots, Eqs. (|2,9( 12.10]) . The simplest nontrivial quantity of that type should be the 
variance var[N n (x\, X2)] of the number of real roots in the vicinity of t = 1. 

Another natural object to study is the variation of the number of real polynomial roots against a 
small change of the vector of real coefficients c = cq, c n -\. More specifically, change c — > c v = c+vh 
assuming the components of the vector b to be i.i.d. standard Gaussian. The parameter v > is used 
here to control the magnitude of the perturbation. By using f v (t) to denote the perturbed polynomial 



an interesting question arises, that is can we generalize the correlation functions Eq. (j3.22j) to the 
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Figure 3: Suppression of the density profile by increasing scaled mean value a from to 8 for a degree 
100 polynomial. 



following parametric correlation functions: 

N^i)(t 1 ,t + A 1 )...N^i)(t l ,t l + A l ) 



Kvi,...,vi(h, -,ti) = lim E 



lim 



A lv ..,A(->0 |Ai|...|A 



(3.23) 



which reflects the change in the positions of the real roots. Similar objects are of interest in the theory 
of random matrices and disordered systems (see |27| and references therein). 

As a preliminary step of our research we evaluated the simplest nontrivial parametric correlation 
function, following the paper [5] and found in the 'local regime': 

K , v {t,t) 1 . 1 

- 1 H — arcsm ■ 



[p(tp v ^y+^ 

It will be interesting to attack the problem in full generality, both for 'global' and especially for 'local' 
regime, where the results are expected to be universal. 

In fact, our initial interest in the properties of random polynomials was stimulated by the fact that 
closely related methods can be applied to study the properties of irregular eigenfunctions in 'quantum 
billiards.' The eigenfunctions ^f(x,y) at the point with the coordinate vector r = (x,y) are solutions 
of the Helmholtz equation: — A^ / (f) = E^(f) ; r € where Q, is a connected compact domain 
with the boundary dQ, and A is the Laplacian. 

Recently Smilansky and collaborators [2S] suggested looking for the eigenfunctions at the point r 
in the following representation: 

L 

V(r,e)= aiMk\r\)e iW (3.24) 

l=-L 

where Ji(x) stands for the Bessel functions, r, 6 are polar coordinates of the observation point, and 
the integer L is given in terms of the wavenumber k and the perimeter D of the billiard boundary 90 
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as L = | [^jp] ■ The complex coefficients ai satisfying a[ = (— Vfa-i are taken to be i.i.d. complex 
Gaussian variables with unit variance, in accordance to The Berry's conjecture. Based on such a 
formula the authors of (2^] managed to calculate the mean number N and variance var{N) of the 
intersection of nodal lines with the billiard boundary d£l. For the Dirichlet boundary conditions along 
the boundary curve dD parameterised as R(0), with < 9 < 2ir, the problem turned out to be 
equivalent to counting the number of real zeros of the function 

L 

u{6)= aiJl(kR(e))e iW 
l=~L 

in the interval 6 € [0, 2ii). Those developments provide an interesting possibility for applying ideas 
and methods from the theory of random polynomials to describe chaotic eigenf unctions. 

Similar methods can be hopefully used to study sensitivity of the nodal lines of eigenfunctions of 
'quantum billiards' with respect to the perturbation of the billiard parameters. For example, one may 
wish to study parametric variations of the quantities involved, with the role of the external parameter 
played by a slight random variation of the boundary curve R{9) or by any other tunable physical 
parameter. 
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